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ABSTRACT 

Aims. In a previous paper we proposed a new model for the emission by amorphous astronomical dust grains, based on solid-state 
physics. The model uses a description of the disordered charge distribution (DCD) combined with the presence of two-level systems 
(TLS) defects in the amorphous solid composing the grains. The goal of this paper is to compare this new model to astronomical 
observations of different Galactic environments in the far-infrared/submillimeter, in order to derive a set of canonical model parameters 
to be used as a Galactic reference to be compared to in future Galactic and extragalactic studies. 

Methods. We compare the TLS model with existing astronomical data. We consider the average emission spectrum at high latitudes in 
our Galaxy as measured with FIRAS and WMAP, as well as the emission from Galactic compact sources observed with the Archeops 
balloon experiment, for which an inverse relationship between the dust temperature and the emissivity spectral index (/8) has been 
proven. 

Results. We show that, unlike models previously proposed that often invoke two dust components at different temperatures, the 
TLS model successfully reproduces both the shape of the Galactic spectral energy distribution and its evolution with temperature 
as observed in the Archeops data. The best TLS model parameters indicate a charge coherence length of ^ 13 nm and other model 
parameters in broad agreement with expectations from laboratory studies of dust analogs. We conclude that the millimeter excess 
emission, which is often attributed to the presence of very cold dust in the diffuse ISM, is very likely caused solely by TLS emission 
in disordered amorphous dust grains. We discuss the implications of the new model, in terms of mass determinations from millimeter 
continuum observations and the expected variations in the emissivity spectral index with wavelength and dust temperature. The 
implications for analyzing of the Herschel and Planck satellite data are discussed. 



1. Introduction 

Early astronomical observations in the far-infrared (FIR) and 
submillimeter (submm) range have revealed the spectral energy 
distribution (SED) of dust emission towards the diffuse interstel- 
lar medium (ISM) to be consistent with a modified black-body 
function of the form 



field, and as a consequence is related to the FIR and submm op- 
tical properties of the dust. Considering a constant gas/dust ratio 
Nh/Mj (with M4 the mass column density), and assuming a sin- 
gle dust temperature along the line of sight (LOS), the specific 
intensity can be written as 



Iy(A)^eeBy(A,Td)xNH 



(1) 



where ly is the specific intensity or sky brightness (energy flux 
per unit wavelength and solid angle). By the Planck function, A 
the wavelength, Td the dust equilibrium temperature, and Nh the 
hydrogen column density. Those measurements were found to 
be consistent with a power-law emissivity function: 



Iv(A) 



with 



4/lo) h- BAA,Td)xMd, 



k(Ao) = eiAo) 



Nh_ 

Md 



(3) 



(4) 



Ee = e(/lo) I — 

\Ao 



(2) 



where eiAo) is the emissivity at the reference wavelength /Iq, and 
jS is the emissivity spectral index, usually taken as/3 - 2 (the so- 
called quadratic law). The quality of the available data generally 
did not allow checking for variations in /3, which was therefore 
assumed to be constant. 

The emission originates in dust radiating at thermal equi- 
librium (at the temperature Td) with the suiTounding radiation 
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where k(Ao) is the dust mass opacity at the reference wavelength 

Aq. 

In solid state physics, within the framework of classical mod- 
els, a constant and temperature-independent emissivity spec- 
tral index value jS = 2 is justified by the asymptotic behavior 
of the absorption cross-section of solids at long wavelengths, 
far from the absorption resonances, such as the silicate absorp- 
tion features at 10-20 fim. However, it was recognized early 
that the quadratic law may not hold for amorphous materials in 
the FIR/submm, since such materials are likely to e xhibit low- 
energy transitions that happen in the FIR/submm (see lAndriessel 
1974). Although big grains (BGs) are known to be amorphous 
in nature in the diffuse ISM (Kemper et al., 2004), no need was 
identified to explicitly incorporate details about their internal 
structure in astrophysical models due to the lack of submm 
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data. However, recent observations have proven that the actual 
FIR/submm dust SED is significantly more complicated than de- 
scribed by Eq.|2l 

Several recent studies have shown evidence of variations 
in the dust emissivity law with wavelength. That the overall 
Galactic SED is flatter than predicted b y the quadratic law 
was clearly shown in the EZRAS data by iReach et aP (Il995h . 
jParadis et al. (2009) studied the spectral dependence of dust 
emission in a collection of molecular clouds and their atomic 
surroundings over the spectral range covered by the DIRBE, 
Archeops, and WMAP experiments and showed that the ob- 
served emissivity law exhibits a break around 500 /im, with a 
steeper index (J3=2A) in the EIR and a flatter index (J3= 1 .5) in the 
submm. More recently a similar trend ha s been apparent in the 
Planck data of the Taurus molecular cloud ( iPlanck coUaborationl 
12011 g|) . Similar submm excess emission h as also been observe d 
towards external low-metallicity galaxies Galliano et al.l (l2005l) . 
as w ell as towards the Small Ma gellanic Cloud in the Planck 
data (iPlanck collaborationll201 lah . The comparison between the 
spectral indices derived for the Galactic solar neighborhood 
(J3 = 1.8), the LMC (P ^ 1.5), and the SM C (/? = 1.2) in the 
Planck data by IPlanck coUaborationl (l2011bllah indicates a gen- 
eral flattening that could be linked to the varying metallicity 
among these three environments. It therefore appears that the 
dust SED is usually flatter than predicted by the quadratic law 
and that the associated emissivity index/? varies with wavelength 
and from object to object, possibly tracing evolution of the dust 
properties with metallicity. 

In addition to the above, the spectral index of the dust emis- 
sion as derived from a modified black-body fit to the data follows 
a systematic trend with dust temperature, the emissivity spec- 
trum appearing flatter for regions with warmer dust temperatures 
and steeper for colder regions. 

- iDupac et alJ (l2003h show variations in p in the EIR in vari- 
ous regions of the ISM, using the PRONAOS balloon-borne 
experiment data. They found variations in the index in the 
range 2.4 to 0.8 for Td between 1 1 and 80 K. 

- iDesert et alj (l2008l) highlighted variations in p in the submm 
towards compact sources in Archeops balloon-borne experi- 
ment data. They found variations in the index in the range 4 
to 1 for te mperatures between 7 and 27 K. 

- IVeneziani et al. (2010) show variations in p between 5 and 

1 in the temperature range 7-20 K by analyzing 8 high 
Galactic latitudes clouds using BOOMERanG observations, 
combined wi t h the IRAS, DIRBE, and WMAP data. 

- iParadis et alJ (|2010) found evidence of T^-P variations using 
the combination of the new Herschel and IRAS data, in the 
inner regions of the Galactic plane. The yS range is 2.6-1.8 for 
temperatures between 14 and 23 K. 

The origin of the flatness of the dust emission spectrum ob- 
served with FIRAS has been explored bv lReach et alJ d 19951) and 
iFinkbeiner et alJ (Il999t) . The latter group propose a phenomeno- 
logical description of the FIRAS data that invokes the coexis- 
tence of a warm (T^ ^ 16 A") and a very cold dust (Td ^ 9 K) 
component (hereafter EDS model). i Reach et aP (11995') deter- 
mined that the FIR and the submm emission were tightly corre- 
lated and therefore did not favor the very cold dust component. 
Moreover, the two component description does not repro- 
duce the observed variations in the emissivity spectral index with 
dust temperature. Indeed the spectral index associated to each of 
the components is set to a fixed value in the EDS model. For 
instance, in their best-fit model, p for the cold silicate-like com- 
ponent is set to 1.7, which does not allow SEDs to be produced 



with a larger apparent emissivity index in the submm/mm do- 
main (see Section lSTl i. 

The two-level systems (TLS ) model (iMenv et all l2007l 
hereafter Paper I) has been developed to describe the FIR to mm 
continuum interstellar dust emission, taking the effect of the dis- 
ordered internal structure of the amorphous dust grains into ac- 
count. It is based on the solid-state physics model developed to 
interpret specific properties of the amorphous solids identified in 
laboratory data. As a consequence, it is expected to apply with 
a high degree of universality, and not to be sensitive to the ex- 
act chemical nature of the dust. The disordered charge distribu- 
tion (DCD) part of the model describes the interaction between 
the electromagnetic wave and acoustic osc illations in the dis- 
ordered charge of the amorphous material dVinogradovl 119601 : 
ISchlomamul 19641) . This charge-disorder is observed on nanome- 
ter scale and is described here by a single charge correlation 
length. The TLS part of the model describes the interaction of the 
electromagnetic wave with a simple distribution of asymmetric 
double-well potential (ADWP), which results from disorder on 
atomic scale: a ADWP can be associated with two close config- 
urations of atoms or group of a toms in the d isordered structure 
dPhilfipsL 119721 Il987t lAnderson et all [T972I) . where the corre- 
sponding atoms have two equilibrium positions and can transit 
from one to the other. Both D CD and TLS phenomena have been 
first applied bv lBoschI (11978) to explain the observed tempera- 
ture dependence of the absorption of some silica-based glasses 
in the FIR/mm. 

In section 2, we present the data we consider for this anal- 
ysis. In section 3, we briefly describe the model in terms of 
the disordered structure, as well as the main parameters of the 
model. In section 4, we explain the methods we used to compare 
the model to astrophysical data and derive standard parameters 
for the Milky Way (MW). In section 5, we present model pre- 
dictions for the standard parameters derived above, such as the 
evolution of the dust emission with temperature, the emissiv- 
ity spectral index as a function of wavelength and temperature, 
and the predicted emissivity in Herschel and Planck photomet- 
ric bands. Sections 6 and 7 are devoted to the discussion and 
conclusions, respectively. 

2. Astronomical data 

In the following we describe the dataset we used in this study. 

2.1. FIR/mm SED of the diffuse Milky Way emission 

The shape of the dust SED in the diflPuse medium is an impor- 
tant constraint on the physics of interstellar dust. In principle, we 
would like to compare the new model with the SED of the most 
nearby and the most diffuse ISM, in order to avoid including 
emission from dense regions where dust properties may differ 
from those in the diffuse solar neighborhood. Also, it is impor- 
tant to analyze regions with minimum temperature mixing along 
the LOS, which may make the interpretations difficult. However, 
waiting for the Planck data to become publicly available, the data 
presently available prevent us from deriving such an SED up to 
millimeter wavelengths, because the low dust emissivity in this 
range and the low dust column density available combine into 
a poor signal to noise ratio (S/N). For instance, th e solar neigh- 
borhood SED presented bv lBoulanger et all d 19961) did not allow 
measuring dust emission in the solar neighborhood for A > l.l 
mm. Instead, we used the emission from the whole MW, but ex- 
clude the central parts of the Galactic disk (\b\ < 6°) in order 
to avoid regions with extensive dust property mixing along the 
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LOS, as well as regions with I240 < 18 MJy/sr, to keep a reason- 
able S/N. Figure [T] shows the region of the sky seen by FIRAS 
at 240 fim that we considered in this work, which accounts for 
13.7% of the entire sky. 

We constructed an average MW SED from 100 yum to 13 mm 
(23 GHz) by merging data from the FIRAS instrument onboard 
the COBB satellite and the WMAP data. We used the full range 
of FIRAS wavelengths (0.1-10 mm or 3000-30 GHz) and the 
5 WMAP bands at 13 mm (23 GHz), 9.1 mm (33 GHz), 7.3 
mm (41 GHz), 4.9 mm (61 GHz) and 3.2 mm (94 GHz). The 
angular resolution of the FIRAS data is assumed to be 7°. The 
angular resolution of the WMAP data is taken from the WMAP 
explanatory supplement to range from 13' to 53' with increasing 
wavelengths. 

For practical purposes, the FIRAS data were first reprojected 
into the all-sky pixelisation scheme HEALPix (Hierarchical 
Equal Area isoLatitude Pixelization) so that it can easily be 
combined with the WMAP data. This was done using the 
method c urrently used for al l Planck ancillary data and de- 
scribed in iParadis et al.l (I2OIII) . The FIRAS data used here are 
available in the HEALPix scheme on the WMAP LAMBDA 
(Legacy Archive for Microwave Background Data) website'. 
The WMAP data were degraded to the FIRAS angular resolution 
(7°) through the convolution with a Gaussian kernel of the appro- 
priate FWHM using the standard convolution tools in HEALPix, 
before extracting the SED. An offset subtraction was done by 
removing an average SED computed over regions of the sky 
with low emission set as I240 < 3 MJy/sr, corresponding to high 
Galactic latitudes. The method used removes any arbitrary addi- 
tive offset in the data that could be set by calibration uncertainty 
or cosmological sources such as the cosmic infrared background 
(CIB). Uncertainties on the SED are computed as the quadratic 
sum of three different uncertainties for each instrument: calibra- 
tion, white noise per pixel, and background. We assumed the 
official calibration and white noise uncertainties given in the ex- 
planatory supplements. Calibration uncertainties are <10% for 
FIRAS, depe nding on the frequen cy range, and less than 1% for 
WMAP data ( Jarosik et al.Ll20I ll) . The background uncertainties 
on the FIRAS and WMAP data have been derived by computing 
the standard deviation over fifty circular regions with a 2° radius, 
taken randomly in the high latitude sky. This uncertainty largely 
predominates over the other ones. 

The obtained SED is shown in Figure |2] and the SED val- 
ues and uncertainties are given in Table[T] The corresponding 
average HI column density is 2x10^' H/cm^. It can be seen 
from Figure 12] that the emission peaks near A = 160 /im, in 
agreeme nt with the FI RAS SED derived at high latitudes by 
IBoulang er et al.l (Il996h . who found an average dust temperature 
of 17.5 K assuming a constant spectral index set to 2. The recent 
analysis of the Planck data shows that a more appropriate value 
for the FIR spectral index in the solar neighborhood is /3-l.E 
dPlanck collaborationl l20IIal) . and using this value over nine 
fields in th e galactic halo leads to a s imilar temperature value 
of 17.9 K (Planc k coUaborationi l20IIfl) and 17.7 K for regions 
with \b\ > 10°. Above A =^ 5 mm, the ISM emission starts in- 
creasing with wavelengths. This contribution does not originate 
in thermal dust emission, but is most probably from a combi- 
nation of synchrotron or free-free emission and/or to anomalous 
microwave emiss ion, which is ofte n attributed to spinning dust 
(see fo r instance iFinkbeiner et aI.L I2004' Planck coUaboratioru. 

l2oTT3dh 



http://lambda.gsfc.nasa.gov/ 



2.2. FIR/submm SEDs of Galactic compact sources 

Before the Planck data become publicly available, the FIRAS 
data are the only data allowing the dust equilibrium temperature 
and the FIR-submm spectral index of dust emission to be mea- 
sures simultaneously over large regions of the sky. This is mainly 
because this requires both FIR data (which measure temperature 
best) and submm data in the Rayleigh-Jeans regime, which con- 
strains the submm emissivity index best. The FIRAS data accu- 
rately sample the peak of the BG emission. A good sampling of 
data in the range 100- 500 fim is crucial for constraining the dust 
temperature, as well as the correlation length, which is a free pa- 
rameter of the TLS model (see Section l372b . Unfortunately, at the 
angular resolution of the FIRAS instrument, the sky distribution 
provides very little handle on dust temperature variations. At the 
FIR AS resolution, the dust temperature varies only from 16 to 21 
K (Reac h et al.Lll995h . assuming a constant spectral index equal 
to 2. Therefore, testing the temperature dependence of the dust 
spectral index requires additional information in the submm over 
wide dust temperature range. 

Balloon experiments have been able to record FIR/submm 
data at a significantly better angular resolution than for the 
COBE and FI RAS data. Both the PRONAOS and th e Archeops 
data analysis (iDupac et all 1200 Uloisert et aUl2008l) show a de- 
crease in the spectral index with temperature. The former study 
used data of the diffuse emission mapped over limited regions 
going from cirrus clouds to dense regions. The latter study was 
carried using observations of compact sources in the MW. The 
variations sample the range from /3 ^ 4 for compact sources 
with dust temperatures as low as 7 K down to y6 ^ 1 for the 
highest temperatures of about 27 K. The interpretation of the ob- 
served variatio ns in the emi ssivity index with temperature is still 
being debated. IShettv et al] (12009 ) argue that the observed anti- 
correlation could result entirely from noise in the data or from 
temperature mixing along the LOS. However, the recent analy- 
sis of the Herschel data in the Galactic plane by P aradis et al.l 
(1201 oh demonstrated that the observed variations could not be 
entirely due to data noise. Similarly, the analysis of the Planck 
data towards cold cores showed that the expected LOS temper- 
ature mixing could not explain the whole observed trend. It is 
therefore likely that the observed steepening of the emissivity 
index with decreasing temperature reveals changes in the opti- 
cal properties of astronomical dust with temperature. Laboratory 
experiments (Boudet et al., 2005; Couceaud et al., 2011) have 
also shown evidence of a similar behavior for analogs of the 
amorphous material likely to compose astrophysical dust grains, 
which tend to favor actual variations in the dust properties with 
temperature. 

Here, we do not intend to fit any given T^-fi relationship. 
Instead, we minimize against the overall Archeops compact 
source spectra, which do not impose following a predetermined 
Td-yS relationship. Archeops observed about 30% of the sky at 
a 12' angular resol ution, in four photometric bands: 550, 850, 
1382 and 2098 um (iBenoit et al.Ll2003 ). The Archeops compact 
sources catalog (iDesert et al.L l2008h contains 302 SEDs, which 
have been fitted using a modified black body. Only 153 fits have 
given a satisfactory ;^'- (using ax^ goo dness-of-fit criteria at the 
2-cr level) when using free Tj and p (iDesert et al.L l2008h . We 
therefore restricted our analysis to these 153 spectra. Since the 
sampling of the compact sources spectra is not evenly distributed 
in dust temperatures, we averaged them in ten bins of tempera- 
ture, between 7.3 and 25.7 K. In this way, we ensure a similar 
weight for the ten spectra. The summary of the binning parame- 
ters, including the bin's central temperature, the average fluxes. 
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Ai^im) 


104.0 


106.5 


109.2 


111.9 


114.9 


117.9 


121.2 


124.6 


128.2 


132.1 


136.1 


140.5 


145.1 


Brightness (MJy/sr) 


17.08 


24.79 


26.28 


30.92 


30.73 


33.27 


32.93 


34.55 


35.72 


36.40 


37.87 


38.86 


39.34 


1-0- (MJy/sr) 


49.38 


15.95 


11.24 


11.66 


7.22 


7.66 


6.48 


6.49 


4.58 


4.37 


3.13 


2.75 


2.28 


A(fim) 


150.0 


155.3 


161.0 


167.1 


173.6 


180.8 


188.5 


196.9 


206.1 


216.2 


227.3 


239.7 


253.5 


Brightness (MJy/sr) 


40.15 


39.90 


41.82 


40.76 


40.86 


40.29 


39.33 


38.13 


36.66 


34.68 


32.51 


30.18 


27.57 


1-0- (MJy/sr) 


2.75 


2.20 


1.39 


1.29 


0.99 


1.50 


1.16 


0.91 


0.53 


0.44 


0.35 


0.28 


0.29 


AQim) 


268.9 


286.4 


306.3 


329.1 


355.7 


386.9 


424.1 


469.2 


525.1 


596.0 


698.1 


816.8 


1002.4 


Brightness (MJy/sr) 


24.92 


22.12 


19.35 


16.62 


14.01 


11.53 


9.37 


7.22 


5.21 


3.66 


2.46 


1.49 


0.79 


1-0- (MJy/sr) 


0.25 


0.35 


0.22 


0.24 


0.24 


0.24 


0.29 


1.01 


0.33 


0.17 


0.12 


0.11 


0.09 


AQim) 


1297.2 




1837.7 


3150.3 


3200.0 


4900.0 


7300.0 




9100.0 


13000.0 


Brightness (MJy/sr) 


0.35 




3.13 


-1.6x10-2 


1.9x10-2 


9.2x10-3 


9.8x10 


-3 


1.1x10-2 


1.5x 


10-2 


1-0- (MJy/sr) 


0.10 




3.10 


0.12 




1.1x10-3 


5.6x10-* 


2.7x10 


-4 


2.3x10-" 


1.8x10-* 



Table 1. FIRAS/WMAP SED of the diffuse MW emission. 



the number of sources, as well as their Archeops names, are pro- 
vided in Tabled 

The uncertainty on the data is the quadratic combination of 
the relative and the absolute uncertainties (15, 17, 21, and 14 % 
from 550 jum to 2. 1 mm). 

3. The TLS model 

In the following, we first remind the reader of the relation be- 
tween dust emission and absorption and then briefly summarize 
the different processes of the TLS model. The full theoretical de- 
scription of the model is given in Paper I. 

3.1. Relation between dust emission and absorption 

The Mie theory is based on the exact resolution of the Maxwell's 
equation for the absorption and scattering of a plane electromag- 
netic wave by a homogeneous spherical particle whose mate- 
rial is characterized by its complex dielectric constant eioj) - 
e'{cL>) + ie"{Lo), with co - IncjA. Solid state physical models of- 
ten analyze the optical properties by the use of the complex sus- 
ceptibility ;^'o('^) = ;k''(<y)-t-iy"(tt>), which in the case of isotropic 
dielectric materials is related to the complex dielectric constant 
by the Clausius-Mossotti equation: 



eicS) - 1 + 



AnxoiMi) 



1- 



Axxai<jj) ' 



(5) 



For spherical grains of radius a and density p, the Mie theory 
gives the following expression for the dust mass opacity: 



k(u)) 



9w le{cj)-l 
— /m ; — - 

pc \ eioj) + 2 



(6) 



The optical properties of bulk materials are also often charac- 
terized by the absorption coefficient a, which is defined as the 
optical depth t in the bulk material per unit length. The absorp- 
tion coefficient a is related to the complex dielectric constant e 
by 



a{u) — 2Im 



'Vi 



(7) 



In general, the complex dielectric constant is wavelength- 
dependent, and as a consequence a{a>) is not proportional to 
K(a)), the ratio K/a is wavelength-dependent. 

However, for dielectric dust materials in the FIR/submm 
range, the value of the imaginary part e" of the complex dielec- 
tric constant e(a>) is negligible with respect to the real part e'(w) 



(details can be found in Paper I). This condition e"((y) <s e'(w) 
leads to 



CO e" 
a(aj) ~ — , 

and therefore: 
K(io) 9 yf? 



(8) 



(9) 



a(w) p(e' -I- 2)2 ■ 

The real and imaginary parts of the complex dielectric con- 
stant are coupled through the Kramers-Kronig relations. It im- 
plies that e"(a>) <& e'(a>) and that the real part of e(a>) can be 
considered constant {de' jdcj - 0). As a consequence, the ratio 
K{(jLi)la{<jS) is almost constant in the FIR/submm range, and the 
slopes of k{cl)) and a{u)) are the same. 

3.2. A disordered cliarge distribution 

The disordered charge distribution (DCD) description of an 
amorphous material shows that in the FIR/mm wavelength 
range, the absorption coefficient due to acoustic lattice oscilla- 
tions is temperature-independent and can be written: 



odcd 



(e + 2Y 1 a 



3 vf m cy/e 



1 



t '^'^ 



(10) 



Here, wo - v,llc with Vt the transverse sound velocity in the 
material and Ic the charge correlation length, c is the speed of 
light, < q^/m > the mean squared charge deviation per atomic 
mass, and e the dielectric constant. It must be pointed out that the 
DCD process presents two asymptotic behaviors for the absorp- 
tion and spectral index, on both sides of wo-The spectral index 
varies between 2 in the short wavelength range {to > wq) and 4 
in the long wavelength range (a> < coo). e, Vt, and < q^/m > is 
likely to vary with material composition and with disorder at sec- 
ond order. In the absence of a clear identification of the material 
composing astronomical grains, we adopt for modeling purpose 
the numerical values used by Bosch ( 1978). The important effect 
for our study is the wavelength profile of the DCD absorption. 
As a consequence only the numerical value of wq is relevant. The 
value of Ic will be derived from wq - v,llc taken the reference 
sound velocity equal to 3x10^ cm.s-' into account. 



3.3. A two-level system 

Originally, the TLS description is based on the assumption that 
the amorphous characteristics of the solid enables close config- 
urations of atoms or groups of atoms, which can be described 
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Bin Tfi, 
(K) 



^100 

(Jy) 



^550 

(Jy) 



(Jy) 



■f'l382 

(Jy) 



^209S 

(Jy) 



nb. compact sources/bin 



Arch, name 



1 8.2 25.1+3.4 



10.1 34.2+4.6 



434.8+66.0 



315.8+47.8 



86.5+15.0 15.1+3.2 3.1+0.6 



65.6+11.3 15.3+3.2 3.0+0.5 



3 


11.9 


99.6+13.4 


380.9+58.5 


81.9+14.4 


21.9+4.6 


4.7+0.8 


4 


13.8 


368.1+49.7 


416.8+63.3 


106.3+18.3 


22.0+4.7 


5.5+0.9 


5 


15.6 


700.8+94.6 


410.1+62.8 


96.2+17.3 


23.2+4.9 


5.9+1.0 


6 


17.4 


1960.7±264.7 


573.3±88.0 


145.6+25.3 


29.8+6.3 


8.1+1.3 


7 


19.3 


3311.4+447.0 


579.6+89.1 


133.0+23.3 


28.6+6.1 


5.9+1.2 


8 


21.1 


1877.5+253.5 


367.3+61.2 


95.5+17.9 


29.9+6.5 


5.6+1.3 


9 


23.0 
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19.3+ 3.3 


10 
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15,41,48,63,64,66,73, 110,124, 132, 

143, 151, 158, 188, 191, 200, 222, 231, 232, 239, 

245, 247, 256, 259, 266, 272, 289, 291, 297 

1, 45, 49, 51, 54, 69, 79, 87, 90, 99, 

105, 106, 111, 137, 138, 139, 141, 145, 146, 150, 

169, 171, 174, 175, 176, 180, 181, 186, 192, 197, 

203, 204, 210, 227, 228, 240, 249, 252, 255, 267, 

270, 274, 276, 284, 295, 296, 298, 

17, 37, 47, 58, 62, 65, 84, 93, 116, 127, 

162, 168, 196, 206, 207, 211, 220, 224, 242, 253 

38, 59, 61, 76, 89, 97, 98, 109, 130, 155, 

157, 160, 166, 173, 185, 189, 216, 229, 251, 260, 

292 

21, 35, 95, 102, 113, 118, 163, 164, 235, 237, 

275,301, 

11, 14, 112, 133, 152, 269, 277, 290, 300 

78, 92, 100, 226, 246 

5, 117, 179,268 

234 

57, 114 



Table 2. Description of the Archeops spectra bins used in this study. 




Fig. 1. Map of the FIRAS 240 yum data used in this work {\b\ < 6° and I240 < 18 MJy/sr). Units are MJy/sr 



by an ADWP. The TLS model describes the interaction between 
electromagnetic waves and the split ground states (which rep- 
resent the two-level systems) of these ADWP that are optically 
active; i.e. they are characterized by a non-zero dipole moment. 
Phenomenologically, because of the amorphous nature of the 
material, each ADWP is local and specific, but in a macroscopic 
analysis, the disorder in the amorphous material can be studied 
in terms of an ordered distribution of TLS with a simple den- 
sity of states. Three interaction mechanisms, which all depend 
on temperature, can occur in such a population of TLS sites: a 
resonant absorption, and two relaxation processes identified as 
"tunneling relaxation" and "hopping relaxation". As opposed to 
resonant phenomena, the electromagnetic field does not cause 
the transition between the two energy levels for relaxation mech- 
anisms. It means that all relaxation processes, regardless of their 
splitting energy Tioj, can participate in the absorption of the elec- 
tromagnetic wave of a given frequence wq, unlike the resonant 
effect where only the TLS characterized by the energy Tiuq can 
absorb. The various effects can be described as follows: 



Resonant absorption dHubbard et al.Ll2003h : 
The resonant absorption of a photon of energy tio) concerns 
only those TLS characterized by the splitting energy tLCj. 
Such a mechanism gives an absorption that decreases with 
the increasing temperature and vanishes when the tempera- 
ture is high enough that the two-levels become equally pop- 
ulated. The resonant absorption can be written as: 



47r2 {e + iy 
3^ 9~ 



a)G(u) tanh(ha)/2kTd), 



(11) 



where k is the Boltzmann constant, e the dielectric con- 
stant, and G{a>) the optical density of state (ODOS). For 
simpli city, and to co mpare our results with a reference, we 
follow iBoschI (119781) and consider a constant ODOS, i.e. 
G(oj) = Go, with Go equal to 1.4x10"^ (COS units). 

Tunneling relaxation: 

In terms of absorption, this process occurs when the phonon 

energy is too low to allow direct transitions above the bar- 



D. Paradis et al.: Far-infrared to millimeter astrophysical dust emission 



rier height between the two levels. In this case, additional 
tunneling is required to perform the transition. The tunnel- 
ing absorption increases with temperature. We derive the fol- 
lowing accurate and relatively simplistic expression for the 
phonon-assisted tunneling absorption, mor e convenient for 
comput ation than the formula obtained by iFitzgerald et al.l 
(I2001allbh : 



^phon 



Go (e + 2y 



3c Vf 



0) F2{(JL>, Td), 



(12) 
(13) 



2kTaJo L T 



Ti E LodrdE 

— sech^(— — ) 3 ^ , 



where ti is the relaxation time defined by ti - 
a E^^ tanh(E /2kTd)- A reasonable val ue for a is 4.2 x 10"^^ 
erg^ s (in the case of soda-silica glass: lB5schL 1 1 9781) . 

Hopping relaxation: 

Hopping relaxation occurs when the phonon energy allows 
a direct jump above the barrier. It gives an absorption that 
increases with temperature: 



^hop — 



--^Go(cA + InTd) dVPiV)- 

9 Jo 1 + w-T 



3c Vf 



^) 



where t is the relaxation time defined by t - tq exp (V/Td)- 
For amorphous silica, tq ^ 1 x 10 '^s (lBoschl[l978l) . ca, an 
additional parameter of tunneling states is given by 



kg 
Ca = In — - + ln4 — I + 



In arctan(x) dx 



(15) 



where A™'" is the minimal energy splitting. Only a lower 
limit of Ca with a value close to 5.8 has been determined 
from laboratory experiments. P(V) is the distribution of the 
TLS barrier height which can be approximated by 



P{V) = CvPviV) with 



Cv = 



1 



VoV^ 



and 



PviV) 



{ exp(-(y - y„)Vv2) for V > v„, 

\ for y < y„ 



(16) 



(17) 



(18) 



such as y„„„ = 50 K, V„ = 550 K and y,) = 410 K. 

Al l physical parameters are set to the values from iBoschI 
(1197 8V and are provided in Table [3] The amplitude of the TLS 
processes with respect to the DCD part is controled by a rela- 
tive intensity parameter denoted A (see Section 13741 1. The total 
absorption {a,o,) deduced from the TLS model is the sum of all 
the processes quoted in this section: 

atot = a DCD + A X (a,-es + aphon + a hop)- (19) 



3.4. Model parameters 

The relevant parameter for the DCD absorption is the charge cor- 
relation length Ic, which controls the wavelength where the in- 
flection point between the two asymptotic behaviors with P - 2 
andyS = 4 occurs. The relative amplitude between the DCD and 
the TLS efifects is not strongly constrained by the model prior to 
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Fig. 2. FIRAS/WMAP diff'use MW spectrum fitted using the 
TLS model: total in red (solid), DCD process in light blue (long 
dash), and TLS processes: resonant absorption in green (dash), 
hopping relaxation in dark blue (dot), and tunneling relaxation 
in orange (dash-dot), with the best-fit parameters (top panel) 
and standard parameters (bottom panel). The gray line (dash- 
dot-dot-dash) corresponds to free-free emission, and the brown 
line (thick dash) to a A" power law. 



the comparison with astrophysical data. The factor that controls 
the amplitude of the DCD absorption depends on the character- 
istics of the material, whereas the amplitude of the TLS effects 
(resonant, tunneling, and partly hopping) is related to the TLS 
density of states. Therefore, to simply characterize the relative 
importance of the TLS and the DCD processes, we use a single 
free parameter of relative intensity denoted A. This means that 
all TLS effects are multiplied by the same intensity A (see Eq. 
[T9] l. The respective contributions of the resonant and of the two 
relaxation processes are set according to theory. The resonant 
effect and tunneling relaxation correspond to a quantum tunnel- 
ing eff'ect. However, the hopping relaxation is a more classical 
effect, since there is no tunneling involved. The characteristics 
of the potential wells to be taken into account in computing the 
absorption in case of hopping or tunneling are different, depend- 
ing on the classical or quantum tunneling process. Therefore, 
we want to keep the possibility of uncoupling these effects, by 
means of the ca coefficient. The dust temperature can be derived 
from the shape of the dust SED in the FIR and also controls the 
amplitude of the TLS relaxation effects. It is therefore a crucial 
parameter of the model, linking the FIR emission to the submm 
behavior of the optical properties of dust. 

Finally, we adopt the following four free parameters for the 
TLS model in this study: r^, Z^, A, and ca. Note that these pa- 
rameters allow only to characterize the shape of the emission. 
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(cm/s) 


< q-/m >' 
erg.cm.g-' 


e*t 


Go 
(CGS units) 


a'- 
(ergls) 


(s) 


^0 

(K) 


y m 

(K) 


mm 

(K) 


3x10" 


8748 


2.6 


1.4x10-' 


4.2x10-'^ 


10-" 


410 


550 


50 



Table 3. Physical parameters used in the TLS model. 
* defined in Section[32l 
' defined in Section [33] 



One additional parameter is required to set the overall amplitude 
of dust emission. Classically this is done by assuming a dust/gas 
ratio and some dust emissivity at a given reference frequency. 
Also, the number of free parameters used here is roughly the 
same as the one needed to achieve a description of the emission 
shape with the FDS model {a\,a2,T2, qilqi, and /i). 



4. Comparison with astrophysical data 

4.1. Determination of tiie best-fit parameters 

The diffuse medium sampled by the FIRASAVMAP SED and 
the compact sources detected in the Archeops data could exhibit 
rather different dust properties. In a first step, we therefore model 
the two datasets separately, in order to illustrate possible differ- 
ences between the two environments. However, as explained in 
Section 12.21 only the Archeops data sample the range of tem- 
perature needed to constrain the TLS model. We also therefore 
perform a combined fit where all data are used, in order to de- 
rive standard parameter values allowing both diffuse medium 
and compact sources to be reproduced. In each case, we pro- 
ceed in the same way, i.e., we use the idl "mpfit" function to 
perform a x^ minimization between the observational data and 
the TLS model. Results of the minimization are given in Table 
m The prediction of the model is color-corrected to match the 
observation in the large photometric filters of each instrument. 



4.1 .1 . Modeling of the FIRAS/WMAP MW spectrum 

The model has been compared with the FIR/mm SED of the dif- 
fuse MW emission. For the longest wavelengths, dust emission 
is not the main contributor, in particular in the WMAP bands. We 
take free-free emission into acco unt using t he free -free template 
at 30 GHz constructed by Di ckinson et al.l ( |2003h . derived from 
the Ha map, and extrapolated in frequency using a power law 
for the specific intensity with ly oc A^-^. In addition to the free- 
free emission, we have incorporated a power law emission that 
can describe synchrotron emission, probably mixed with other 
types of emission such as spinning dust. The slope of this power 
law (a) is left as a free parameter, since contamination in this 
wavelength range is not well known. The brightness level of the 
free-free combined with the power law emission has been nor- 
malized to the 23 GHz data. Figure|2](top panel) shows the result 
of the modeling. From Table|4]we can see that the dust tempera- 
ture is robustl y derived, and the v alue of 17.5 K is similar to the 
results of Boulanger et alJ (Il996h . who analyzed the emission in 
the solar neighborhood. However the other parameters are not 
well constrained using only the FIRAS/WMAP SED. The A pa- 
rameter controls the intensity of all TLS effects, whereas the ca 
parameter allows us to increase/decrease the emission coming 
from the hopping relaxation. To reproduce only the mm flatten- 
ing in the spectrum deduced from the FIRAS data, the model re- 
quires an increase in the emission in the submm domain, which 
can be explained by an increase of either A or ca. These two 



parameters are therefore highly degenerate in this fit configu- 
ration. The 1-cr uncertainty on the correlation length indicates 
that using only the FIRAS/WMAP data is not sufficient for con- 
straining this parameter, which controls the slope of the emission 
spectrum. Indeed, changing Ic is equivalent to modifying the ap- 
parent emissivity spectral index/?, and/3 is known to be naturally 
antic orrelated with dust temperature. Therefore, the only way 
to constrain this parameter is to compare the model to spectra 
at various temperatures, which means that the modeling of the 
Archeops spectra is important to derive Ic (see Section |4.1.2| |. 
The specific intensity is best fitted by a power law of the form 
ly oc A" with a - 1.45. Variations in the synchrotron index have 
been found in our Galaxy with values ranging from 0.6 in star- 
forming complexes to 1.2 in the interarm regions. Our value is 
higher than the observed a values for synchrotron alone, which 
could indicate that this emission is not only due to synchrotron 
emission, but may also include a significant contribution from 
spinning dust emission. Although spinning dust emission has 
not been shown at high la titude, some studies indicate that it is 
widespread in the Galaxy (" Miville-Deschenes & Lapachell2005l : 
Planck collaboratioa.201 Id) . 

4.1 .2. Modeling of the Archeops spectra 

In this section, we perform the minimization on the Archeops 
spectra, averaged in ten temperature bins. The IRAS 100 jum 
data were also included in the spectra. Since the wavelength do- 
main is limited to A ^2 mm, no free-free or spinning dust emis- 
sion were included in the fits. For each set of parameters [l^ A, 
Ca], the model is computed for temperatures between 5 K and 
35 K, with a step of 0.5 K, and adjusted to the spectra for the 
best dust temperature for each spectrum, by looking for the low- 
est ;t'^ value. Therefore, we determine the best set of parameters 
that is able to reproduce all Archeops spectra by only changing 
the dust temperature. Results of the minimization are presented 
in Figure [3] and Tabled A comparison between the temper ature 
derived from a modified black-body fit jDesert et al.L 120081) and 
the temperature deduced from the TLS model are also given in 
Figure [3] Whereas a modified black-body fit samples the range 
8.2-24.8 K, the TLS model shows temperatures from 10 to 17.5 
K. This discrepancy is essentially due to the mm excess which 
can bias the temperature determination when using a modified 
black-body fit. In Section 1531 we discuss the temperature esti- 
mates depending on the method used. 

The parameters are better constrained than in Section 14.1.1! 
with lower 1-cr uncertainties, especially for l^ with a value of 
5. 11 ±0.09 nm. In Figure [3] we can see that the flattening of 
the spectra with increasing dust temperature (from top to bot- 
tom panels, and from left to right) is mostly consistent with the 
hopping relaxation process. However, the spectra with the low- 
est dust temperature indicate that the Archeops data at longer 
wavelengths and at 2 mm in particular, show a resonant absorp- 
tion on the same order of magnitude or higher than implied by 
the hopping relaxation process. Also, we note that the tunneling 
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Best fit parameters: lc=5.11nm; A=3.86; cA=1333 
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Fig. 3. Spectra of Archeops compact sources (diamonds), binned in temperature, fitted using the TLS model, with the best-fit pa- 
rameters (this page) and the standard parameters (next page). The contribution of the different processes is given in the caption of 
Figure |2] Squares correspond to the model after color correction. The average temperature per bin deduced from a single modified 
black-body fit (from Desert et al.,2008), and the temperature deduced from the TLS model are given in the top right of each panel. 
The description of the Archeops spectra is given in Section 2.2 and in Table 2. 



relaxation contribution is much lower than the other processes. 
The minimization of Galactic compact sources favors a lower in- 
tensity of the TLS resonant process which results in the increase 
in the amplitude of the hopping relaxation, as compared to the 
modeling of the FIRAS/WMAP spectrum. In other words, the 
decrease in the A parameter requires a higher value of c^. 



4.1 .3. Modeling of the FIRAS/WMAP and Archeops spectra 

To derive general properties of Galactic dust, we performed a 
combined fit of the diffuse medium and compact sources SEDs. 
Best-fit parameters, denoted "standard" parameters are given in 
Table |4] Figures |2] and |3] also show the results of the model- 
ing when using the standard parameters. Since the dust temper- 
ature is only derived from the FIRAS/WMAP SED, we would 
expect a similar value as in the combined fit. However, because 
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Standard parameters: lc=13.40 nm; A=5.81; cA=475 
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Fig. 3. Continued. 



the FIR spectral shape of the BG emission also depends on /c, 
changes in Ic produce changes in Td, since the two parameters 
are degenerate. We note a decrease in /(. with a value of 13.4 nm, 
compared to the value derived using only the diffuse medium 
emission (/^ = 23.05 nm). This explains that the best temper- 
ature estimate (Jd - 17.26 K) does not fully match the value 
obtained when fitting only the FIRAS/WMAP data {Td = 17.53 
K). Actually, the Ic, A, and ca parameter values are intermedi- 
ate, compared to the values obtained when fitting independently 
the diffuse medium emission on one side and compact sources 
on the other. We note that using the same set of TLS parameters 
for all Archeops compact sources does not reproduce the long 



wavelength data well in the case of the warmest spectra as deter- 
mined by a modified black-body fit (7/^,7=23.0-24.8 K.). In that 
case, the TLS model underestimates the brightness at 1382 and 
2098 iim. 

In Figure |2] we do not see any significant difference in the 
modeling of the FIRAS/WMAP spectrum when using the best- 
fit or standard parameters. However, when looking at Figure 
[3] the use of the standard parameters has some non-negligible 
effects on the fit modeling. For instance, the fitting of the 
coldest spectrum shows an overestimate of the brightness at 
long wavelengths, which was not the case with the best-fit pa- 
rameters. Also, the warmest spectra are described less well 
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Table 4. Parameters of the TLS model. 





Td(K) 


4 (nm) 


A 


c^ 


a 


reduced y^ 


Galactic diffuse emission (FIRAS/WMAP) 


Parameters 
l-cr 


17.53 
±0.02 


23.05 
± 22.70 


9.38 
± 1.38 


242 
± 123 


1.45 
±0.15 


1.95 


Galactic compact sources emission (Archeops) 


Parameters 
l-cr 


— 


5.11 
±0.09 


3.86 
±0.13 


1333 
±68 


_ 


1.45 


Galactic diffuse and compact sources emission 


Standard parameters 
1-0- 


17.26 
±0.02 


13.40 
± 1.49 


5.81 
±0.09 


475 
±20 


1.55 
±0.11 


2.53 



when adopting the standard parameters. This could indicate that 
dust properties are different in the diffuse medium observed 
with FIRAS/WMAP and in compact sources observed with 
Archeops. 



4.2. Discussion of fit resuits 



Parameters giving the best match between the TLS model and 
the Galactic emission spectra correspond to a dust temperature 
ofTd =17.26 K, a correlation length close to /c = 1 3 nm, an inten- 
sity of the TLS processes with respect to the DCD one of about 
A =6, and an intensity parameter for the hopping relaxation of 
ca =475. The dust temperature characterizing the diffuse inter- 
stellar medium i s quite close to the value of 17.5 K derived by 
iBoulanger et alj (|T996) at high Galactic latitudes. We note that 
the correlation length is within the same order of magnitude as 
the minimal size of the smallest grains in the BG population. The 
largest grains could result from aggregation of grains with a size 
close to the correlation length. The value of the A parameter is 
reasonable considering the paucity of qu antitative experimental 
data in this field (essentiallv lBoschll 19781) . In solid-state physics, 
the phenomenology of these effects is different, and these effects 
could result from distinct types of disorder in the material. The 
moderate value found for the parameter A indicates that the rela- 
tive intensity of the DCD and TLS pro cesses has the same order 
of magnitude as in the case studied bv lBoschI (Il978h . 

For the physical properties related to defects in the grain 
material, the measured values of defect concentration and in- 
duced defects can easily vary over several orders of magnitude. 
Therefore the derived A value suggests a similar origin for the 
DCD and TLS processes, and therefore justifies the use of a 
similar density of state. The value of c^ is significantly higher 
than the minimum value of 5.8 deduced from laboratory ex- 
periments, at low temperature. Once again, in a material with 
a known chemical composition, the magnitude of the physical 
effects caused by defects in the material are likely to vary over 
several decades. The value derived from observational data is in 
that sense fully acceptable. This value indicates that the hopping 
relaxation is more efficient than the two other TLS processes, 
related to tunneling and therefore highly dependent on the dis- 
tance between the two potential minima of the TLS. However, 
we note that each TLS process is characterized by its own emis- 
sion profile, dependent on the parameters for each site. For in- 
stance, modifying V^, Vq, and Vmin would induce variations in 
the emission profile of the hopping relaxation effect. As a con- 
sequence the Ca value would have to be modified. 



5. Model predictions 

5.1. Spectrai siiape evoiution witti temperature - Emissivity 
spectral index 

In Figure |4] we present the predicted model SED evolution us- 
ing the standard parameters taken from Table|4]for different dust 
temperatures. The resonant absorption and the hopping relax- 
ation process have different behaviors with temperature: the hop- 
ping relaxation increases with increasing temperature, whereas 
the resonant absorption decreases. The tunneling relaxation can 
be neglected over the range of temperature explored here, since 
it appears to be systematically dominated by one of the other 
TLS processes, and does not induce any specific variations in 
the dust emission spectrum. The millimeter excess emission, as 
well as the flattening of the spectrum with temperature, is mainly 
produced by the hopping relaxation. As a consequence, the ex- 
cess increases with temperature, and the transition between the 
DCD and TLS processes shifts towards shorter wavelengths with 
increasing temperature since the intensity of the hopping relax- 
ation increases. 

One of the most important aspects of the TLS model to 
be explored is the variation in /3 with temperature. We com- 
pared the TLS to the FDS model number 7 (see Table 2 in 
Finkbein er et al.L Il999h . The same method has been applied to 
each model. We generated a set of emission spectra for temper- 
atures between 10 and 35 K corresponding to temperatures de- 
rived using the TLS model (hereafter Ttls ) in one case, and to 
the warm component (T2) of the FDS model, which is related 
to that of the cold component (see Eq. 14 in Finkbeine r et all 
1999), in the other Each spectrum has been integrated in the 
IRAS 100 jum and Archeops photometric bands to produce sim- 
ulated SEDs. We then fitted these fake SEDs at different tem- 
peratures by a single modified black body (Td-yS fit), assum- 
ing the same absolute uncertainties as for the actual IRAS and 
Archeops data (see Section l272l l. Results of the T^-fS fits are pre- 
sented in Figure |5] The TLS model clearly highlights an inverse 
relationship between the spectral index and the temperature of 
the fit (here after Tfit) in agreement with the dust behavior in 
iDesert et alj ([2008), whereas the FDS model shows the opposite 
behavior, i.e. an increase m/3 with temperature. We also note that 
the FDS model is not able to predict y6 values higher than 1.6 in 
the wavelength range observed with Archeops. Using standard 
parameters, the TLS model predicts /3 variation between 2.3 and 
1 in the wavelength domain going from 100 jum to 2 mm, for Tfit 
between 10 and 70 K. 

Models such as the FDS model, which do not exhibit any 
intrinsic variation in the apparent spectral index with tempera- 
ture may still explain the obervations in principle, provided an 
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adequate temperature mixing along th LOS or within the beam 
sufficiently alters the emitted SED. 

First, we note that, since temperature mixing can only 
broaden the emitted SED, the intrinsic yS should be larger than 
the range of observed ones. Sinc e many of the Archeops poin t 
sources and Planck cold cores (iPlanck coUaborationi 1201 let) 
have apparent/? larger than 3-4, this would require a model with 
a very steep intrinsic (3 value. Apart from the TLS model pro- 
posed here, which can produce this steep behavior under specific 
circumstances, there has been no proposal for dust models that 
allows such high yS values. 

Second, the effect of temper ature mixing has been analyzed 
in several studies. . Malinen et al.l(l2011 ). for instance, have inves- 
tigated this effect through 3D MHD simulation, taking both full 
radiative transfer and a dust model with no intrinsic f5 variations 
into account. They considered both cases, with and without em- 
bedded stars heating the surrounding ISM, assumin g reasonable 
luminosity and spatial distribution for the sources. . Paradis et al.l 
(|2009|) studied the emiss ivity variations i nduced by an ISRF 
strength mixture using the lDale et al.l (1200 ih concept of a power- 
law distribution of dust mass subjected to a given heating inten- 
sity (dMdiXjsRF ": Xj^j^p-dXjsRF)- The a parameter has been 
tuned in the range a=l (active star-forming regions) up to a-2.5 
(diffuse medium). Both methods give similar results. 

For cases with no embedded source, temperature mixing 
can reproduce neither strong wavelength-dependent variations 
in the emissivity nor the spread of /5 values observed. Indeed, 
iMalinen et al.l (l201lh show that the model corresponding to that 
case actually produces a weak correlation between Tj and /?, 
which is the opposite of the anticorrelation that is observed. 
In cases of star-fo rming regions, when discrete stellar heating 
is included in the iMalinen et al.l (1201 Ih model, most LOS in 
their simulation show a similar correlation trend. However, those 
LOS towards discre te objects (which co rresponds to Dale's a- 
parameter ^L25 in iParadis et al.L 120091) present a wavelength- 
dependent variation in emissivity, leading to an apparent Td-y6 
anticorrelation with a shape similar to what is observed toward 
some regions. Those map pixels affected by temperature-mixing 
around heating sources, ho wever, have a ppare nt temperatures 
higher than 14 K following Malinen et al.l (1201 ih and, as akeady 
pointed out, /? values always lower than the ones characterizing 
the SEDs without temperature mixing. 

As shown in Table |2] most of the Archeops point sources 
considered here have temperatures lower than 14 K, down 
to 8 K, and high spectral index values. This shows that they 
are most likely not affected by temperature mixing and are 
actually closer to co l d cores , as detected in the Planck data 
(Planck coUaborat ioni 1201 lei) . It is there observed that such 
sources exhibit narrow SEDs. For a single dust temperature 
along the LOS, their narrowness is only compatible with 
spectral index values higher than 3. And a strong distribution of 
temperature would require still higher spectra index values to 
fit the observed SEDs. For the above reasons, we do not favor 
a model without intrinsic /3 variations and temperature mixing 
along the LOS to simultaneously explain the diffuse MW SED 
and the Archeops point sources, and propose that the TLS model 
with the parameters derived here could reasonably explain both. 

To extend the analysis on the /3 variations, we present in 
Figure|6]how the predicted ^0 changes as a function of both tem- 
perature and wavelength, computed as dlnee/dlnA. For temper- 
atures below 15 K, the spectral index is expected to exceed 2 
and reach 2.6 at Td=10 K, around 2 mm. Then/? decreases with 
increasing wavelength and gets as low as 0.5 at longer wave- 




100 



1000 

X {\im) 



10000 



Fig. 4. Prediction of the FIR/mm emission at different temper- 
atures, deduced from the TLS model using standard parameters 
(see TableHJl. The contribution of the different processes is given 
in the caption of Figure |2] 



lengths, where dust emission is usually dominated by other emis- 
sions such as free-free and synchrotron. Temperatures higher 
than 20 K and 100 K lead to /3 values less than 2 and 1.7 at 
A > IQQ fim, respectively. We note that at 2 mm, in the range 
10-40 K, the spectral index roughly varies from 2.6 to 1 . To sum- 
marize, the yS values predicted by the TLS model vary strongly 
both with wavelength and dust temperature. 

In Figure |7] we show the evolution of /? in three ranges of 
wavelengths ([100 /zm-2 mm], [100-550 yum] and [550 jjxn-l 
mm]) for different couples of the parameters [A, Ic]. The ca value 
is fixed to the standard value of 475. For a given yS and tempera- 
ture, one can estimate the correlation length of the dust grains, as 
well as the intensity of the TLS effects. Regardless of the value 
of the parameters. Figure |7] highlights a decrease m/3 with tem- 
perature going from 10 to 55 K for most [A, Ic] couples. The 
/3 variations are more pronounced as A increases. The behavior 
of /3 with temperature does not change for Ic > 30 nm, regard- 
less of the A value. In that case, the DCD part has reached the 
asymptotic behavior in A^^, and only the TLS effects can in- 
duce a decrease in /3 with Td- Low values of Ic associated with 
low values of A produce the highest /? values, which can reach 
;0=3.4-3.5 at low temperatures in the three wavelength ranges. 
The DCD asymptotic behavior in /I "* has not been reached in 
Figure |7] since the TLS effects never vanish completely and Ic 
values around 0. 1 nm would be required, which is not realistic 
since the interatomic distance is ^0.2-0.3 nm range. High values 
of A induce low values of /3 at high temperatures. The value of 
/3 can be as low as 0.55 between 550 /im and 2 mm, for temper- 
atures higher than 25 K, in case of A = 10 and Ic = 1 nm. The 
same values of A and Ic, however, reproduce yS ^ 1.4, for Td-55 
K in the range 100-550 ;um. 
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(Finkbeiner et al. 1999) 
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Fig. 5. Prediction of the emissivity spectral index as a function 
of temperature, deduced from a modified black-body fit on emis- 
sion spectra at various temperatures, in the IRAS 100 //m and 
Archeops bands. The emission spectra have been generated us- 
ing the TLS (black) (using standard parameters, see Table|4|i and 
FDS (gray) models. 
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Fig. 6. Prediction of the emissivity spectral index deduced from 
the TLS model using standard parameters (see Table |4|, as a 
function of wavelength and temperature, with overplotted con- 
tours. 



Fig. 7. Prediction of the emissivity spectral index in the range 
[100 A/m-2mm] (top panel), [100-550 //m] (middle panel), [550 
fiin-2 mm] (bottom panel), deduced from the TLS model as a 
function of dust temperature, for different parameter couples [A, 
Ic (nm)] assuming Ca = 475. 



5.2. Dust emissivity 

One of the most important aspects of a given FIR/mm dust emis- 
sion model is to predict the dust emissivity and its variations with 
wavelength and temperature. This point is pivotal both for infer- 
ring gas masses and for separating foreground emission from the 
cosmic microwave background (CMB). Here, we quantify how 
adopting the predictions of the TLS model affects the dust mass 
determination obtained from FIR and submm measurements. 
The reference emissivit y value is 1x10"^^ cm^ /H at 250 jum in 



the solar neighborhood (Boulan ger et al.L [19961) . This value has 
been derived for a dust temperature of 17.5 K. Besides varia- 
tions with wavelength, the TLS model predicts variations with 
temperature. We impose the model to recover the Galactic emis- 



sivity reference at 250 iim from iBoulanger et al.l (Il996h . for a 
temperature of 17.5 K. Emissivity predictions have been inte- 
grated for different temperatures in the IRAS, Herschel PACS 
and SPIRE, Planck HFl and LFl channels. The color correction 
has then been applied to each instrument. Values of the predicted 
emissivities are given in Table |5] Figure [8] shows the emissivity 
predictions for temperatures between 5 K and 100 K, overlaid 
with a A^^ power-law emissivity. We note that between 5 and 17 
K, the emissivity variations follow a A^^ power law in the range 
100 /im-1 mm. At higher temperatures, the emissivity spectra 
become flatter than a modified black-body with y6=2, starting in 
the FIR domain and flattening even more with wavelength. We 
observe that emissivities at 5 and 10 K only differ for A > 4 
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perature derived using the TLS model and the Td-/3 fit is modeled 
using a polynomial function, such as 



TtLS ^ZjJnTji,. 



(20) 



Fig. 8. Predicted emissivity with the TLS model using the stan- 
dard parameters given in Tabled as a function of wavelength for 
different temperatures. The dashed line shows a A^^ emissivity 
power law for comparison. 



mm. For Td >10 K, we can see that the difference between the 
curves increases significantly with wavelength and temperature. 
Whereas the discrepancy between predictions at 5 K and 100 K 
is about a factor 2 at 250 /vm, it is about a factor 35 at 2 mm, and 
^100 at 4 mm. 

Therefore, when assuming no dependency of the emissiv- 
ity spectrum with either wavelength or temperature, one can 
significantly bias the mass estimate, especially when derived 
from submm/mm data and extrapolated to shorter wavelengths. 
Indeed, when assuming an emissivity decrease in A^^ in the 
submm/mm domain, one would overestimate the dust temper- 
ature and underestimate the emissivity in the FIR, and therefore 
would overestimate the dust mass. For instance, this is the case 
for those SMC molecular clouds whose masses have been de- 
rived from mm data and found t o be systematica lly higher (twice 
on average) than virial masses (iBot et al.Ll2007h . 

5.3. Temperature estimate 

Dust temperature is often derived using a modified black-body 
fit with Td and p taken as free parameters, or with p set to a 
fixed value. The latter temperature does not necessarily equal 
the physical temperature obtained using the TLS model, as al- 
ready discussed in Section lSTI In particular, the temperature de- 
duced from a Td-y6 fit is highly dependent on the wavelength 
range considered. To compare the fitted temperature with the 
physical one, we first generated several emission spectra at dif- 
ferent temperatures between 10 and 60 K using the TLS model. 
We then integrated the TLS model predictions in the IRAS 100 
yum band (crucial for estimating the dust temperature), Herschel 
PACS 160 //m, and SPIRE 250, 350, and 500 fim bands. We 
then fitted each spectrum using a least-square fit method. We 
also performed the analysis by integrating the TLS model in 
the IRAS 100 /im channel, and Planck HFI 350, 550, and 
850 fim bands (857, 545, and 353 GHz), before fitting the 
spectrum with a modified black-body model. For both cases 
(IRAS -Herschel bands and IRAS-Planck bands) we have also 
made the test using a Td-j6 fit with fi - 2. For each band, we 
took the calibration uncertainties into account, i.e. 13.5% for 
IRAS (a bsolute uncertainty of the new generation of IRAS maps 
[Miville- Deschenes & Lagac he. 2005), 20% for PACS, 15% for 
SPIRE, and 7% for HFI. The relationship between the dust tem- 



The best polynomial coefficients derived from a fit are given in 
Table|6] 

The comparison between the temperature estimates is shown 
in Figure |9] The behavior varies significantly when deriving the 
temperature in the 100-500 jum range, or in the 100-850 /im 
range, and whether a free (3 is used. Up to 25 K, we observe 
good agreement between Ttls and Tot, no matter which model 
is used, whereas it is not the case for higher temperatures. When 
using a fixed (free) p, Tflt is systematically lower (larger) than 
Ttls ■ Also, the discrepancy between the temperatures is greater 
when enlarging the wavelength domain to derive Tot, and when 
P is left as a free parameter. This result highlights that the longest 
wavelengths {A > 550 jim) should not be used to determine the 
dust temperature from a single modified black-body fit, since 
they are affected by significant submm/mm emission excess in 
the TLS model. In addition, one should be careful when using a 
Td-yS fit, especially when yS is a free parameter of the fit, and for 
Tfit > 25 K. In that case, Equationl20land Table |6] should be used 
to derive the actual physical dust temperature from a temperature 
deduced from a Td-/? fit. 

We note that if the temperature estimate is performed on a 
wavelength range that extends up to 2 mm, as is the case when 
using the Archeops data for instance, the discrepancy between 
Tfit and Tjls is shifted, and roughly starts and increases for tem- 
peratures over 18 K, instead of 25 K (see temperatures presented 
in Figure [3]l. 

6. Discussion 

One important aspect of dust emission is to be able to subtract 
the foreground emission. According to the TLS model, emission 
from large ISM grains made of amorphous material is signif- 
icantly more complicated than predicted by a single modified 
black-body emission with a single emissivity index p. Actually, 
the model predicts emissivity index variations as a function of 
both temperature and wavelength. Therefore, within the frame- 
work of this model, it is important to take into account both 
dependences for component separation. Extrapolation of emis- 
sion from the FIR to mm wavelengths with a single modified 
black body with a constant yS could produce a wrong estimate of 
dust emission in the CMB wavelength domain. The above vari- 
ations are strictly related to the amorphous nature of the mate- 
rial composing the grain, and should not be mistaken with varia- 
tions that can result, for instance, from modifications of the grain 
structure, such as ch a nges in the grain siz e or dust aggregation 
(IStepniketaLl 120031: iParadis et aU. l2009l) . which can occur in 
dense regions of the ISM. 

However, we caution the reader that the interpretation of the 
model prediction may become inaccurate when high tempera- 
tures (^100 K) are concerned. Indeed, in the current version of 
the TLS model, only the fundamental energy of the TLS states is 
considered, and the model requires improvemen t by adding the 
effect of excited states in the two-level systems dGromov et al.L 
120111) . 

The TLS model predicts a wavelength and temperature de- 
pendence of the spectral index that characterizes the emission 
and absorption properties of the amorphous dust in the FIR and 
mm range. Coming from solid states physics, it should apply 
with a high degree of universality to all amorphous materials. 
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Table 5. Predicted emissivities in the IRAS, Herschel, and Planck photometric bands, deduced from the TLS model with the standard parameters given in Table |4l as a function 
of temperature. The emissivities have been normalized to 1 x 10"^^ cm^/H, at 250 fim for a temperature of 17.5 K. Note: Table 5 is also available in machine-readable form in 
the electronic edition of A&A. 
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Table 6. Polynomial coefficients (see Eq. 
modified black-body fit. 



relating the dust temperature derived from the TLS model and those derived from a 




Fig. 9. Temperature derived using the TLS model (Ttls ) using 
the standard parameters given in Table |4] as a function of tem- 
perature derived from a modified black-body fit with a free /3 
pai-ameter (Tfit) in the IRAS 100 fj.m - Herschel 160, 250, 350, 
and 500 /urn channels (black stars), and in the IRAS 100 jum - 
Planck 350, 550, and 850 j-im channels (gray stars). The trian- 
gles show the temperatures derived from a modified black-body 
fit with 13-2, using the IRAS-Herschel channels (black) and 
the IRAS-Planck channels (gray). The dashed line corresponds 
to Tfit = Ttls • The continuous lines shows the best polynomial 
fits. 



Obviously, the exact profile of the emission spectrum depends 
on parameters that define the disorder, i.e. the amorphous struc- 
ture of a material. Indeed, the amorphous state is not well defined 
as crystals, and atoms can be organized in a variety of local ar- 
rangements giving rise to a diversity of structures on short and 
medium length scales. But the emission profiles present some 
common characteristics: in the low-temperature range a spec- 
tral index value of the order of 2 in the FIR range, followed by 
variations in the yS towards higher or lower values at longer wave- 
lengths. As the temperature increases, the absorption is enhanced 
and is characterized by /? values that decrease. These specificities 
can be observed in Figure|6] for the given sets of parameters that 
fit the observational data. However, this general behavior is also 
observed in laboratory experiments performed on amorphous 
materials and grains of astrophysical interest (lAgladze et aU 
Il996t iJVIennella et ail [1998; 'Boudet et al.', '2005'). In particular, 
recent studies on three analogs of amorphous silicate Mg2Si04, 
IVlgSiOa, and CaMgSi206 reveal that this temperature and wave- 
length dependence of the absorption is observed on all the sam- 
ples, but disappe ar when the same sa mples are annealed until 
crystallization (Coupe audet al.L 1201 ih . These analyses indicate 
that the Tj-^S variations are likely to result from intrinsic dust 
properties that can be reproduced for the first time by a dust 
emission model such as the TLS one. 



7. Conclusions 

The theoretical TLS model described in a previous work (Paper 
I) and based on the amorphous structure of grains was com- 
pared to astrophysical data, such as the FIRAS/WMAP and 
Archeops data. The FIRAS/WMAP spectrum represents the dif- 
fuse interstellar medium in the stellar neighborhood, whereas the 
Archeops data characterize dust properties in a variety of com- 
pact sources, where a significant inverse relationship between 
the dust temperature and the emissivity spectral index had been 
proven. We performed a x^ minimization to determine standard 
values for the parameters of the TLS model selected to capture 
the spectral and temperature variations in the model. These free 
parameters are the dust temperature {Td), the charge correlation 
length (Ic) that controls the wavelength where the inflection point 
between the two yS = 2 and j3 = 4 ranges occurs, the intensity 
factor of the TLS processes (A) with respect to the DCD effect, 
and the intensity factor of the TLS/hopping process (ca). Results 
indicate that emission in the submm/mm is dominated by the 
hopping relaxation. According to the model, the BG emission in 
the FIR/mm domain depends on wavelength and temperature, 
which is fundamental both for dust mass determination from 
FIR/submm measurements but also for component separation. 

Using the best-fit parameters (Td = 17.26 K, /c=13.4 nm, 
A=5.81, and ca=475) allowing reproduction of both the emis- 
sion from the diff'use medium and the compact sources, the 
model predicts significant/? variations for temperatures between 
5 and 100 K, with a maximum value of 2.6 at 2 mm. The 
TLS model is presently the only model able to predict /? vari- 
ations with temperature and wavelengths, as observed in both 
observational and laboratory data. The dust emissivity can be 
seriously underestimated if its variations with temperature and 
wavelength are not taken properly into account, generally induc- 
ing overestimates of the dust mass. Similarly, the comparison 
between temperatures derived from the TLS model (Ttls ) and 
those derived from a Tj-yS fit (Tfit), between 100 and 500 fixn, or 
100 and 850 fj.m, are in good agreement for temperatures below 
25 K. Above this value, the derived temperatures differ signif- 
icantly from each other, in particular when emission far from 
the emission peak is considered in the fit. Indeed, a modified 
black-body fit is too simplistic to reproduce the mm emission 
excess, which increases with temperature. However, our study 
enables the possibility of deducing Ttls from a given Tgf We 
also predicted dust emissivities in the IRAS 100 /urn, Herschel, 
and Planck bands for Ttls between 5 and 100 K, which are use- 
ful for comparison with the Planck and Herschel data. For in- 
stance, the provided tables and figures can be used to derive the 
charge correlation length and the intensity of the TLS processes 
from the Herschel and Planck data. 

The TLS model is the first physically-motivated model able 
to reproduce both diffuse medium and compact sources, ob- 
served with FIRAS/WMAP and Archeops. In the future, the 
Planck data will be used to systematically compare the model 
prediction with submm data. 
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